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We calculate the quarkonium formation time in relativistic heavy-ion collisions from the space- 
time correlator of heavy quark vector currents in a hydrodynamics background with the initial 
nonequilibrium stage expanding only in the longitudinal direction. Using in-medium quarkonia 
properties determined with the heavy quark potential taken to be the free energy from lattice 
calculations and the fact that quarkonia can only be formed below their dissociation temperatures 
due to color screening, we find that Y(1S), Y(2S), Y(3S), J/ip and ip' are formed, respectively, 
at 1.2, 6.6, 8.8, 5.8, and 11.0 fm/c after the quark pair are produced in central Au+Au collisions 
at the top energy of Relativistic Heavy Ion Collider (RHIC), and these times become shorter in 
semi-central collisions. We further show, as an example, that including the effect of formation time 
enhances appreciably the survivability of T(1S) in the produced hot dense matter. 
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I. INTRODUCTION 


Studying quarkonium production in relativistic heavy- 
ion collisions is useful for verifying the existence of a 
quark-gluon plasma (QGP) in the collision and under¬ 
standing its properties as first suggested by Matsui and 
Satz 1|_. Studies based on various theoretical mod¬ 
els |2lTl0l] have shown that the suppressed production 
of quarkonia observed in relativistic heavy ion collisions 
at SPS, RHIC, and LHC |TD ITbt is mainly due to their 
dissociation by the thermal partons in the produced 
quark-gluon plasma, particularly when its temperature 
is high [E], E2 . Although these studies all give reason¬ 
able descriptions of the experimental data, various as¬ 
sumptions and model parameters were introduced. One 
of these is the time for quarkonia production from ini¬ 
tial hard nucleon-nucleon collisions. This time is relevant 
for determining the survival probability of a quarkonium 
in the hot dense matter as the heavy quark pair before 
forming the quarkonium is not likely to be dissociated in 
a medium either due to the color screening effect or by 
scattering with the thermal partons as shown in a recent 
study by one of the authors based on the color evapo¬ 
ration model fl8| . The survivability of a quarkonium is 
thus low if its formation time is short, as it is produced 
in the hot dense matter, and high if its formation time is 
long, as it is more likely to be formed outside the matter. 
For quarkonium production from a nucleon-nucleon col¬ 
lision in vacuum, its formation time is known to be not 
short compared to the time for the production of a heavy 
quark pair [l9l - l2lj |. The quarkonium formation time in 


heavy ion collisions is, however, not well determined. In 
some studies it is taken to be its value in the vacuum, 
while in other studies it is assumed to be the same as the 
thermalization time used in the hydrodynamic approach. 
It has also been treated simply as a parameter in some 
studies. 

Using the space-time correlator of heavy quark vector 
currents and the heavy quark potentials extracted from 
the lattice QCD, we have shown in a recent study that 
the presence of a QGP makes the quarkonium forma¬ 
tion time longer, and this effect becomes stronger with 
increasing temperature of the QGP HU- This study is, 
however, carried out for a constant QGP temperature, 
thus neglecting the effect due to the rapid decrease in 
temperature. In the present study, we extend our previ¬ 
ous calculations to include this effect by using the time- 
dependent QGP temperature obtained from a hydrody¬ 
namic model. 

This paper is organized as follows: In Sec. IE we briefly 
review the quarkonium formation time in the vacuum and 
in a QGP. We then describe in Sec. HE the quarkonium 
formation time in a hydrodynamic background for rela¬ 
tivistic heavy-ion collisions. Finally, a summary is given 
in Sec. IIVI 


II. FORMATION TIME OF QUARKONUM IN 
VACUUM AND IN QGP 

By using the dispersion relation, the space-time cor¬ 
relator of the heavy quark vector current operator can 
be decomposed into the propagators of physical states 

as Enm 
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where 


D{ Sl r 2 =-x 2 ) = -^. K i(V st) (2) 

is the relativistic causal propagator in the coordinate 
space with K\ being the modified Bessel function and 
r being the Euclidean proper time j23|. At finite tem¬ 
perature, the factor tanh( v / s/2T) should be multiplied to 
the right hand side of Eq . JTJ, but its effect is negligible 
for the present study I2 ll 1261 . 

Each propagator is weighted by the imaginary part 
ImII(s) of the heavy quark polarization function, which 
has contributions from both resonances and the contin¬ 
uum states: 


ImII(s) 


® e< 2 i / /nu2 r ,/2 

(VS-M<) 2 +r?/4 

e 2 

+ ~r-0{\/s — \fhii)- (3) 


In the above, eQ is the electric charge of a heavy quark 
in units of e; |-i/> i (0)| is the coordinate space wavefunction 
of the resonance at the origin; and ^/Sth is the threshold 
energy for the continuum states. The mass distributions 
of resonances are taken to have the Breit-Wigner form 
with a peak at Mi and the width Tj, and the continuum 
states are described by a step function for simplicity. 

Substituting Eq. (0 into Eq. ([]]) leads to following con¬ 
tributions from resonances and continuum states to the 
heavy quark space-time correlator: 


IT(t) 

hicont ) 


q M 2 

-r-mo )? Ki {M z t), 

7T Z T 

Sen r°° 

Ki(x)dx ■ 


(4) 


These expressions show that the contribution from the 
continuum states dominates at early times and is grad¬ 
ually taken over by that from resonances [2ll . |22| . The 
formation time of the ground state is defined as the time 
when its contribution dominates over all other states. Al¬ 
though r is the Euclidean time, the above definition leads 
to the same formation time or length of a quarkonium as 
in other approaches 0, H3] • 

To calculate the formation time, we first consider the 
fraction of the ground state in the correlation function, 
F(t) = IIo(r)/n(r), and then differentiate it with re¬ 
spect to r, P(t) = dF(r)/d,T. The function P(r) indi¬ 
cates how rapidly the ground state becomes dominant 
in the correlation function. The formation time is then 
defined by the expectation value 


(vforrn ) 


f dr tP(t) 
J dr P(r) ■ 


(5) 


For the formation times of excited states, the contribu¬ 
tions from lower-energy states are first subtracted from 
the correlation function. The fraction of the first excited 


state, for example, is given by F(t) = ni(r)/(n(r) — 
IIo(t)) such that the contribution from the first excited 
state, III (t) , eventually dominates [2l], [22j . 

According to Eq. 0, the quarkonium formation time 
is determined by the three parameters, Mj, |-0(O)|, and 
^/sth- For charmonium formation in vacuum, its mass 
is known and the threshold energy can be taken to be 
twice the D meson or B meson mass. The quarkonium 
wavefunction at the origin is related to the dielectron 
decay width of the quarknoium by 


167ra 2 eQ 

Mf 


IV'(O)! 2 . 


( 6 ) 


Using the values determined from the experimental data 
on e + e _ —> QQ for these parameters, the formation times 
of J/tj) and T {IS) in the vacuum are 0.44 and 0.32 fm/c, 
respectively 121 1 | . 

For quarkonium formation in a QGP, the values of 
Mi , 1-0(0) |, and -y/sth can be determined from solving 
the Schrodinger equation: 


2 m Q -V 2 + V (r, T) 

m Q 


Mr,T) 


Miipi(r, T), (7) 


using the heavy quark potential V (r, T) and heavy quark 
mass mQ in the QGP. According to our recent study 
based on the QCD sum rules the free energy poten¬ 
tial from lattice calculations is the appropriate heavy- 
quark potential for quarkonium at finite temperature if 
the quark mass is taken to be its vacuum value. Tak¬ 
ing the continuum threshold energy to be 2 ttjq + V(r = 
oo, T), i.e., the energy of the heavy quark pair when they 
are infinitely apart, we have evaluated the values of Mi, 
|0(O)|, and y/sth as functions of the QGP temperature. 



FIG. 1: (Color online) Formation times of charmonia (left 
panel) and bottomonia (right panel) in a QGP as functions 
of the temperature of the QGP fisS. 122*]. 

Figure [T] shows the formation times of charmonia and 
bottomonia in a QGP as functions of its temperature [H], 
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[22| . We can see that they are longer than their values 
in vacuum and increase with temperature. These results 
are consistent with the picture that it takes a longer time 
for the heavy quark pair to form a quarkoium when its 
radius becomes larger with increasing temperature [22j . 


III. QUARKONIUM FORMATION TIME IN 
RELATIVISTIC HEAVY-ION COLLISIONS 

The above study can be extended to relativistic heavy- 
ion collisions by taking into account the change of the 
temperature of the QGP as the quarkonium forms from 
the heavy quark pair. This is achieved by using the 2+1 
ideal hydrodynamic model, assuming the boost invari¬ 
ance and a zero viscosity for the QGP, which we have 
previously used in studying quarkonia production in rel¬ 
ativistic heavy ion collisions [23, |30| . Specifically, we as¬ 
sume that for heavy ion collisions at the top energy of 
RHIC, hydrodynamics is applicable after an initial ther- 
malization time of about 0.6 fm/c [2(| [3(3]. To apply 
the method described in the previous section to the for¬ 
mation of quarkonia in the nonequilibrium matter before 
the thermalzation time, we assume for simplicity that the 
modification of the quarkonium formation time is inde¬ 
pendent of the momentum anisotropy of the matter [3lj 
and is determined by an effective temperature that is re¬ 
lated to the entropy density of the matter through the 
lattice equation of state [H] that is used in our hydrody¬ 
namic model. Further assuming that the initial nonequi- 
librium matter expands only in the longitudinal direc¬ 
tion, which is supported by studies based on the Color 
Glass Condensate picture in which the pre-equilibrium 
state in relativistic heavy-ion collisions is composed of 
color electromagnetic fields [H, 0], the entropy density 
is then inversely proportional to time. 

Since quarkonium production in heavy ion collisions is 
from initial nucleon-nucleon binary collisions, we use an 
average temperature for the hot medium by weighting the 
local temperature T(x,y,r) with the number of binary 
collisions Aginary ( x i y) i 11 the transverse plane, that is, 

T (j^ = I (,X(, y '/'(•r-//.r)-Vi )illa ,- v (.r.//) ^ 

f dxdy A^bi nar y (x,y) 

We calculate Aginary using the Glauber model with the 
nucleon distribution in a nucleus given by a Woods-Saxon 
function and the empirical nucleon-nucleon scattering 
cross section as described in Ref. [35j |. Shown in Fig. [2] 
by the red dotted lines in upper and lower panels are 
the average temperature of produced matter in Au+Au 
collisions at ^/snn = 200 GeV for the 0-5% and 50-60% 
centralities, respectively. Assuming that the quark pair is 
produced at r = 1/(27 tiq), where toq is the heavy quark 
mass and is taken to be 1.5 and 4.75 GeV for charm and 
bottom respectively, it is seen that the initial tempera¬ 
ture in collisions at the 0-5% centrality is so high that 
only the continuum states can be formed. 



FIG. 2: (Color online) Average temperature T, the fractions 
of bottomonium states F, and their derivatives P as functions 
of time in Au+Au collisions at ^/snn = 200 GeV for the 0-5 
% (upper panel) and 50-60 % (lower panel) centralities. 


Once the temperature drops to the dissociation tem¬ 
perature of Y(1S) state, the resonance begins to ap¬ 
pear near the continuum threshold energy and then com¬ 
pletely separates from it. The dashed and solid lines in 
Fig. [5] are, respectively, the fraction F of each bottomo¬ 
nium state in the correlator for the bottom quark vector 
current and its derivative P with respect to r. The black, 
blue, and green dashed lines, shown from left to right in 
the figure, correspond to the IS, 2S, and 3S bottomonium 
states, respectively. Each resonance starts to appear at 
its dissociation temperature, which is 2.5 T c for the IS 
state, 1.1 T c for the 2S state, and 1.0 T c for the 3S state, 
where T c = 158 MeV is the phase transition temperature, 
when the lattice free energy potential is used. In colli¬ 
sions at the 0-5 % centrality, the system reaches the dis¬ 
sociation temperature at r = 0.42 fm/c for the IS state, 
at 6.4 fm/c for the 2S state, and at 8.4 fm/c for the 3S 
state, and these values change to 0.11 fm/c, 2.4 fm/c, and 
3.4 fm/c, respectively, in collisions at the 50-60 centrality. 
We see that it takes longer for the higher excited state 
to form. We also see that bottomonia are formed later 
in central collisions and also take longer than in semi¬ 
central collisions. The formation times of the IS, 2S, 
and 3S states are, respectively,1.2, 6.6, and 8.8 fm/c in 
Au+Au collisions at the 0-5% centrality and 0.7, 2.6, and 
3.8 fm/c in collisions at the 50-60 % centrality. Although 
the excited states of bottomonium start to form at later 
times, they complete their formation more quickly than 
the ground state. This can be understood as follows. 
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Since the heavy quark and antiquark pair produced in a 
hard collision are initially close in space, it takes time for 
them to separate to a distance that is comparable to the 
size of a quarkonium in a medium. It thus take longer for 
the heavy quark pair to form excited quarkonium states 
due to their larger sizes. However, these excited states do 
not start to form before the system reaches their dissocia¬ 
tion temperatures when the heavy quark pair are already 
sufficiently separated, thus requiring less time for them 
to subsequently form the excited quarkonium states. 



t (fm/c) 


sions is so high that the heavy quark and antiquark pair 
produced in a hard collision are screened from each other, 
there is still a possibility for them to form a bound state 
if the duration of the high density stage is short, because 
they are initially very close to each other in space. How¬ 
ever, if the duration of high energy density is long, it is 
then unlikely that the initial heavy quark and antiquark 
pair can form a quarkonium as they are far apart at the 
formation time, such as for the excited states of bottomo- 
nia and all charmonia in central collisions as well as the 
excited states of both bottomonia and charmonia in semi¬ 
central collisions. In this case, quarkonia production in 
heavy ion collisions will be dominated by the regener¬ 
ation from thermalized heavy quarks and antiquarks in 
the QGP. 



FIG. 3: (Color online) Average temperature T, the fractions 
of charmonium states F, and their derivatives P as functions 
of time in Au+Au collisions at ^/snn = 200 GeV for the 0-5 
% (upper panel) and 50-60 % (lower panel) centralities. 

Figure [3] shows the charmonia formation time in 
Au+Au collisions at ^/snn = 200 GeV for the 0-5% and 
50-60% centralities, respectively. The dashed and solid, 
and dotted lines are, respectively, the fractions F of char¬ 
monia, their derivatives P, and the average temperature 
T in collisions corresponding to the two different central¬ 
ities. Since the dissociation temperature of J/ip is much 
lower than that of T(1S), it appears later at both cen¬ 
tralities. With the lattice free energy as the heavy quark 
potential, the dissociation temperatures of J/ip and ip' 
are about 1.14 and 0.94 T c , respectively. The time for 
this to happen is 5.6 fm/c for J /ip and 10.4 fm/c for ip' in 
collisions at the 0-5 % centrality and is 2.0 fm/c for J/ip 
and 4.4 fm/c for ip' in collisions at the 50-60 centrality. 
Both charmonium states appear after their dissociation 
temperatures and take, respectively, 5.8 and 11.0 fm/c to 
form in collisions at the 0-5% centrality and 2.2 and 5.1 
fm/c in collisions at the 50-60% centrality. 

Although the initial energy density in heavy-ion colli- 


FIG. 4: (Color online) Thermal decay width of T (IS) as a 
function of temperature (left panel) and its survival probabil¬ 
ities from thermal decay with and without including the for¬ 
mation time as functions of collision centrality (right panel). 

After its formation, a quarkonium can be dissociated 
by thermal partons. Its survivability in a heavy ion col¬ 
lision is given by 

O /o“ drP ( r )exp(-/ T Tcr (T')dT') 

^ COO j r>( \ ■ W 

Jo drP(T) 

if we neglect the thermal dissociation below T c as its ef¬ 
fect is much smaller than that above T c . In the above, 
P(t) is the fraction of quarkonium formed per unit time 
at t as defined above Eq. d5j) and r c is the time when the 
quark-gluon plasma stage ends. The dissociation width 
T(r) of the formed quarkonium at time r in the above 
equation depends on the temperature of the quark-gluon 
plasma at that time and its dissociation cross section by 
thermal partons. For the latter, it can be calculated from 
the pQCD with massive thermal partons as in Ref. [ljf| 
and used extensively in Refs. 0, [lcj, 0, HH, 37]. As an 
example, we show in the left panel of Fig. 0] the dissocia¬ 
tion width of Y(1S) as a function of temperature calcu- 
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lated with the lattice free energy as the potential between 
the bottom quark and antiquark pair in the quark-gluon 
plasma. In the right panel of Fig. [1] we show by filled 
squares the survival probabilities of T(1S) in Au+Au col¬ 
lisions at y / SNN = 200 GeV for both centralities of 0-5% 
and 50-60%. Compared to those for T(1S) of zero for¬ 
mation time (filled circles), i.e., the T (IS) is formed as 
soon as the temperature drops to its dissociation temper¬ 
ature, including the formation time increases appreciably 
the survivability of T(1S) in relativistic heavy ion colli¬ 
sions, particularly in central collisions. 

IV. SUMMARY 

The formation time of quarkonium is one of essen¬ 
tial ingredients used in studies on quarkonium produc¬ 
tion in relativistic heavy-ion collisions. To calculate the 
formation time in the quark-gluon plasma produced in 
relativistic heavy-ion collisions, we have extended the 
method which uses the space-time correlator for the 
heavy quark vector current and its dispersion relation 
to the case that the produced quark-gluon plasma is de¬ 
scribed by a 2+1 ideal hydrodynamic model. Before the 
initial thermalization time where hydrodynamics is not 
applicable, we have assumed that the hot dense matter 
expands only in the longitudinal direction and its effect 
on quarkonium formation can be described by an effective 
temperature. We have found that the quarkonium states 
begin to appear as the hot dense matter expands and 
cools to their dissociation temperatures, and the times 
to complete their formations in central Au+Au collisions 
at the top energy of RHIC are at 1.2, 6.6, 8.8, 5.8, and 


11.0 fm/c for T(1S), Y(2S), Y(3S), J/4> and if', respec¬ 
tively. In semi-central or peripheral collisions, both the 
start time of the formation and the formation time of a 
quarkonium are shorter than in central collisions. 

The increased formation time of quarkonia in the hot 
dense matter is expected to affect their survivability in 
relativistic heavy ion collisions. We have found, for exam¬ 
ple, that including the formation time increases apprecia¬ 
bly the survivability of Y(1S). This effect depends, how¬ 
ever, on the magnitude and temperature dependence of 
the thermal decay width of a quarkonium, with a strong 
temperature dependence giving a larger effect and a weak 
temperature dependence giving a small effect. Unfortu¬ 
nately, both the magnitude and temperature dependence 
of the quarkonium thermal decay widths are presently 
not well determined as very different values have been 
used in the phenomenological studies of qaurkonium pro¬ 
duction in relativistic heavy ion collisions. 

Finally, since the formation time of T(IS') in relativis¬ 
tic heavy ion collisions is not much longer than the initial 
thermalization time, a better treatment of the effect of 
the initial nonequilibrium matter than that used in the 
present study is needed. 
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